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'Tn Abstract. We compute the exact relaxation rate of the partially asymmetric 

exclusion process with open boundaries, with boundary rates opposing the preferred 
Q direction of flow in the bulk. This reverse bias introduces a length scale in the system, 

1 '> at which we find a crossover between exponential and algebraic relaxation on the 
Cy coexistence line. Our results follow from a careful analysis of the Bethe ansatz root 
C^ structure. 
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1. Introduction 

The partially asymmetric exclusion process (PASEP) \U2i is one of the most thoroughly 
studied models of non-equilibrium statistical mechanics |3]-[6] . It is a microscopic model 
of a driven system [T] describing the asymmetric diffusion of hard-core particles along 
a one-dimensional chain. In this paper we will consider a finite chain with L sites. At 
late times the PASEP exhibits a relaxation towards a non-equilibrium stationary state. 
In the presence of two boundaries at which particles are injected and extracted with 
given rates, the bulk behaviour at stationarity is strongly dependent on the injection 
and extraction rates. The corresponding stationary phase diagram as well as various 



physical quantities have been determined by exact methods [3j-[6|[8 13 . 

More recently exact dynamical properties such as relaxation rates have been 
derived. On an infinite lattice where random matrix techniques can be used, considerable 



progress has been made, see e.g. 14 17 . For the PASEP on a finite ring, where particle 



number is conserved, relaxation rates were obtained by means of the Bethe ansatz some 



time ago 18 21 . The dynamical phase diagram of the PASEP with open boundaries 



was studied in 22-24 , again using Bethe ansatz methods. 

In the present work we extend some of these results to the case of reverse bias, 
where the boundary parameters strongly oppose the bias in the bulk hopping rates. We 
will in particular focus on the coexistence line where the effects are most pronounced. 

1.1. The partially asymmetric exclusion process 



a 



A A A 



Y 5 

Figure 1. Transition rates for tire partially asymmetric exclusion process. 

We now turn to a description of the dynamical rules defining the PASEP on a one 
dimensional lattice with L sites, see Figure [TJ At any given time t each site is either 
occupied by a particle or empty. The system is then updated as follows. 

Particles attempt to hop one site to the right with rate p and one site to the left 
with rate q. The hop is prohibited if the neighbouring site is occupied, and at sites 1 
and L these rules are modified in the following way. If site i = 1 is empty, a particle may 
enter the system with rate a. If, on the other hand, this site is occupied, the particle 
will leave the system with rate 7 or attempt to hop to the right with rate p. Similarly, at 
i = L particles are injected and extracted with rates S and (3 respectively, and attempt 
to hop to the left with rate q. 



It is customary to associate a Boolean variable Tj with every site, indicating whether 
a particle is present (tj = 1) or not (tj = 0) at site i. Let |0) and |1) denote the 
standard basis vectors in C^. A state of the system at time t is then characterised by 
the probability distribution 

\P{t)) = J2Pir\t)\T), (1.1) 

T 

where 

L 

|T) = |ri,...,ri) = (g)|r,). (1.2) 

1=1 

The time evolution of \P{t)) is governed by the aforementioned rules, which gives rise 
to the master equation 

^|P(t)) = M|P(t)), (1.3) 

where the PASEP transition matrix M consists of two-body interactions only and is 
given by 

L-l 

M = J2 ^^^""^^ ® Mfc ® j(^-*^-i) + mi ® J(^-^) + /(^-^) ® rriL. (1.4) 

fc=i 

Here /*■*''-' is the identity matrix on the A;-fold tensor product of C^ and, for 1 < A; < L — 1, 

Ma. : C2 ® C2 ^ C^ 
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The terms involving mi and ttil describe injection (extraction) of particles with rates a 
and 6 (7 and (3) at sites 1 and L respectively. Their explicit forms are 



^1= ' I ^.= 7 o ■ (1-6) 





The transition matrix M has a unique stationary state corresponding to the 
eigenvalue zero. For positive rates, all other eigenvalues of M have non-positive real 
parts. The late time behaviour of the PASEP is dominated by the eigenstates of M 
with the largest real parts of the corresponding eigenvalues. 

In the next sections we set p = 1 without loss of generality and determine the 
eigenvalue of M with the largest non-zero real part using the Bethe ansatz. The latter 
reduces the problem of determining the spectrum of M to solving a system of coupled 
polynomial equations of degree 3L — 1. Using these equations, the spectrum of M can 
be studied numerically for very large L, and, as we will show, analytic results can be 
obtained in the limit for large L. 



1.2. Reverse bias 

The reverse bias regime corresponds to the boundary parameters opposing the direction 
of flow in the bulk. In this paper we achieve this by considering the regime 

g < 1, a,/3<7,(5. (1.7) 

With g < 1 the bias in the bulk is for diffusion left to right, but setting a, (3 ^ 7, 5 
strongly favours particles entering from the right and exiting on the left. Traditionally 
the reverse bias regime was considered only to exist for the case a = /3 = 0, see e.g. 13|f 



2. Stationary phase diagram 

In this section we briefly review some known facts about the stationary phase diagram. 
It will be convenient to use the following standard combinations of parameters flo] , 

«a,7 = ^ (^a,7 ± y <7 + 4^7) , V^,^ = 1 - q ~ a + -f. (2.1) 

In order to ease notations we will use the following abbreviations, 

'^ — '^a,7) " — ^p,Si ^ — ^0,7; d — ^/3,<5- (2-2) 

The parameters a, 6, c, d are standard notation for the parameters appearing in the 



Askey- Wilson polynomials Pn{x] a, 6, c, d\q) used for the PASEP sationary state 25 . 



2.1. Forward bias 



The phase diagram of the PASEP at stationarity was found by Sandow |T(^ and is 
depicted in Figure |2j The standard phase diagram is understood to depend only on the 



parameters a and h deflned in (2.2) rather than p,q,a, j3,'y,5 separately, see e.g. the 
recent review IB]. However, we will show below that the picture may be more nuanced 
in the regime where a, /3 ^ 7, 6. 

Many quantities of interest, such as density proflles and currents, have been 
calculated for particular limits of the PASEP. For the most general case these have 



been computed by 25 , 26 . 



2.2. Reverse bias 

Much less is known for the case of reverse bias. The stationary state normalisation and 



current ioi a = (3 = have been computed by 13 . It is found that the current decays 
exponentially as 



It was argued in 13 that non-zero values of a and (3 would allow for a right flowing 



current of particles to be sustained, thus destroying the reverse bias phase. However, in 

I Note that the authors of [131 consider q > 1 and 7 = (5 = but we can simply translate their results 
to our notation using left-right symmetry. 
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Figure 2. Stationary state pliase diagram of the PASEP. Tiie liigh and low density 
plrases are separated by the coexistence line (CL). The maximum current phase (MC) 



occurs at small values of the parameters a and b defined in (2.2 1. 



the following we will see that when a and /3 are 0{q^) for some ?7i > 0, the system still 
feels the effects of a reverse bias in its relaxation rates. It would be interesting to know 
whether such effects survive in the stationary state. 



3. Relaxation rates — approach to stationarity 

The slowest relaxation towards the stationary state at asymptotically late times t is 
given by e^^*, where £i is the eigenvalue of the transition matrix M with the largest 
non-zero real part. 

It is well-known that the PASEP can be mapped onto the spin- 1/2 anisotropic 
Heisenberg chain with general (integrable) open boundary conditions flOl 11 . By 



building on recent progress in applying the Bethe ansatz to the diagonalisation of the 



Hamiltonian of the latter problem 27 29 , the Bethe ansatz equations for the PASEP 
with the most general open boundary conditions can be obtained. Recently these 
equations have also been derived directly using a (matrix) coordinate Bethe anstaz 
in the PASEP 30,31 . The Bethe ansatz equations diagonalise the PASEP transition 
matrix. By analysing this set of equations for large, finite L the eigenvalue of the 
transition matrix with the largest non-zero real part has been obtained in the forward 
bias regime 22 24 . In this regime, particles diffuse predominantly from left to right 
and particle injection and extraction occurs mainly at sites 1 and L respectively. We 
will briefly review these results. 



3.1. Bethe ansatz equations 



As was shown in 22-24 , all eigenvalues £ oi M other than £^ = can be expressed in 



terms of the roots Zj of a set of L — 1 non-linear algebraic equations as 
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where 



^o = a + /3 + 7 + (5=(l-g) 



1 — ac 



l-hd 



(l + a)(l + c) (l + 6)(l + d) 



and the "bare energy" e{z) is 

{q-lfz 



eiz 



il-z)iqz-l) 



;i-g) 



z — 1 qz — 1 



The complex roots Zj satisfy the Bethe ansatz equations 



1 — Zi 



K{z,) = n 



qzj - zi q^ZjZi - 1 



^^. Zj - qzi ZjZi - 1 



Here K{z) is defined by 



K{z) 



{z + a){z + c) {z + b){z + d) 

{qaz + l){qcz + 1) {qbz + l){qdz + 1) ' 



and a, b, c, d are given in (2.2). 



L-l. 



(3.1) 



(3.2) 



(3.3) 



(3.4) 



(3.5) 



3.1.1. Analysis in the forward bias regime The analysis of the Bethe equations for the 
forward bias regime proceeds along the following lines. One first defines the counting 
function, 

1 1 ^"^ 

iYLiz) = g{z) + -g^{z) + -J^ ^izi, z), (3.6) 



where 



L^ 



1=1 



9{z) 
9h{z) 



In 



z{l — qzY 



.■w£(i-4!li+ta 



z + a 1 + c/z 
1 + qaz 1 + qcz 



+ ln 



K{w,z) 



1-^2 

z + b 1 + d/z 
1 + qbz 1 + qdz 

1 — qz/w 1 — q^zw 



\\i{w) — In 



1 — qw/ z 1 — wz 



(3.7) 



We note that the choice of branch cuts in this definition of K{w.i z) differs from 24 



but is numerically stable for large L. A further discussion on the choice of branch cuts 



in K is given in Section 4.2 



Using the counting function, the Bethe ansatz equations (3.4) can be cast in 
logarithmic form as 

yL{zj) = '^I,, j = l,...,L-l, (3.8) 



where Ij are integers. Each set of integers {Ij\ j = 1, ... ,L — 1} in (3.8) specifies a 
particular (excited) eigenstate of the transition matrix. 

In l24l, the logarithmic Bethe equations were solved numerically to find the 
structure of the root distribution (an example of such a solution will be seen in the 
next section). By comparing the resulting eigenvalues with brute force diagonalisation 
of the transition matrix, the integers Ij corresponding to the first excited state can be 
found for small values of L. The integers were found to be 

Ij = -L/2+j for j = l,...,L-l. (3.9) 



Assuming (3.9) holds for all L, (3.8) can be turned into an integro-differential equation 
for the counting function Yl{z), which can be solved asymptotically in the limit of large 
L. Exact asymptotic expressions were obtained for the first excited eigenvalue in the 
forward bias regime in all areas of the stationary phase diagram except the maximal 
current phase. We recall here the expression for the coexistence line: 

gi= \~\ _f . +0{L-% (3.10) 

L^ [a ^ — a) 

3.2. Dynamical phase diagram 

The dynamical phase diagram for the PASEP resulting from an analysis in the regime 
g < 1 is shown in Figure [3] It exhibits the same subdivision found from the analysis of 
the current in the stationary state: low and high density phases (a > 1 and 6 > 1) 
with dynamic exponent z = 0, the coexistence line {a = h > 1) with diffusive 
{z = 2) behaviour and the maximum current phase {a,b < 1) with KPZ {z = 3/2) 
universality. Furthermore, the finite-size scaling of the lowest excited state energy of 
the transition matrix suggests the sub-division of both low and high-density phases into 
four regions respectively. These regions are characterized by different functional forms 



of the relaxation rates at asymptotically late times 24 , see also 132] . 



4. Bethe ansatz in the reverse bias regime 

The reverse bias regime corresponds to taking parameters a, b large, and c, d towards 
—1 (note that —l<c,d< 0). This change, and the resulting changes in the structure 
of the root distribution, requires a reformulation of the counting function, which will be 



discussed in Section |4.2| and [Appendix A[ But first we describe the changes to the root 
structure. In the following we address solutions on the coexistence line, a = b, and to 
simplify computations further also set c = d. 

Moving from the forward to the reverse bias regime, the structure of the Bethe 
roots changes. From a forward bias root distribution (Figures |4l 5(a)), increasing — c 




Si 



Figure 3. Dynamical phase diagram of the PASEP in the forward bias regime 
determined by the lowest excitation E\. The horizontal axes are the boundary 
parameters a and h (2.2 1 and the vertical axis is the lowest relaxation rate. The 



latter goes to zero for large systems on the coexistence line (CL) and in the maximum 
current phase (MC). The curves and lines correspond to various crossovers in the low 
and high density phase, across which the relaxation rate changes non-analytically. 



pushes out a complex conjugate pair of roots along with the single root on the positive 



real axis (Figure 5(b) ). As — c increases further, the complex conjugate pair pinch in to 



the real axis until some point where they "collapse" (Figure 5(c) ) - one root returns to 



close the contour on the positive real axis leaving an isolated pair of real roots at 



— c ± e 



-v^L 



v^ >0. 



(4.1) 



In fact, there is a pair of isolated roots for all A; G N such that — g'^c falls outside the 
contour of complex roots. 
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Figure 4. Root distribution in the forward bias regime for L ~ 40, with parameters 
a = 6 = 11, g = 0.3 and c — —0.07. The dashed line marks the position of — c and the 



region in the box is expanded in Figure 5(a) below. 
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Figure 5. Leading part of tiie root distribution for L = 40, a — h — 11, q = 0.3, 
witli c = —0.07 (forward bias), c = —0.088 (intermediate stage), and c = —0.1 (reverse 
bias). Tlie dashed line shows the position of — c. 



For L large, the contour crosses the real axis at |24| 

z*{a) = ^ '- ~ -, 4.2 

la a 

where the approximation is for a large. Therefore, the integer m defined by a and c 

through the inequalities 



g'"c < - < -g'"-'c, 
a 



(4.3) 



is equal to the maximum number of pairs of real roots as L — )■ oo. While we take our 
parameters to have generic values, it would be of interest to analyse the system at the 
resonant values of a, c where the inequalities become equalities, see e.g. p53]. We define 



4x10- 



-4x10-" 



0.2 



0.4 



0.6 



'"% 



-4x10- 



u 



2.x 10-'" 4.x 10-'" 6.x 10-'" 8.X10- 



0.8 



1.x 10 



Figure 6. Complete root distribution (top) and zoomed in on the contour and first 
few pairs of isolated roots (bottom) for L = 100, a = 3.09, c = —0.9, q — 0.3 and 
TO = 20 (at this resolution only five real pairs are visible in the top figure). 



10 



g'"-'", 



a = q- ^a, (4.4) 

then, for example, we will reach a maximum m = 20 pairs of real roots with the 



parameters 

a = 3.09, c=-0.9, g = 0.3. (4.5) 

The corresponding root distribution is displayed in Figure [6j 

We denote by m! the number of pairs of real roots for a particular value of L, with 
m! < m. Then the number of roots on the contour is L — 2m' — 1, which is plotted in 
Figure [T] for the parameters (4.5). There are three regions to consider: 

(i) Isolated roots region: L — 2m' — 1 small and almost all roots occur as isolated pairs; 
(ii) Crossover region: The contour of complex roots begins to form, m' ~ -m; 
(iii) Asymptotic region: The number of roots on the contour, L — 2m' — 1 is large. 
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Figure 7. Number of complex roots in tiie contour (L — 2m' — 1) against L. Only 
beyond L = 40 does the contour of complex roots begin to form. 



4.I. Isolated roots region 

We first consider the case where L — 2m' — 1 is small and assume all roots are real, 
a situation similar to that encountered in the recent calculation of the PASEP current 
fluctuations i34j. The m' pairs of isolated roots are 



-q''c±e 



-utL 



k = 0, 



>m - 1, 



(4.6) 



and the remaing root is close to —(f^c. Substituting these roots into the expression for 

)pp 

;i-g)- 



the eigenvalue (3.1) and dropping exponentially small parts, the sum telescopes, leaving 



^1 



(is) 



g 



m—\ 



+ a 



(1 - q)q 



m' — l 



-qc 



+ 



-q^c 



1 + q'^'c 1 + g™'+ic 



• (4.7) 



The L dependence enters through m' and implies that the eigenvalue decays as q^^"^ to 
a constant, which is 0{q"^). This differs markedly from the asymptotic expression in 



the forward bias regime (3.10). 



11 



4-2. Asymptotic region 

We now analyse the third region, where the number of roots on the contour is large. 
As we are interested in the asymptotic behaviour, we will assume we have reached the 
maximum number of pairs of real roots {m' = m). We label the roots on the contour 

l...L-2m-l, (4.8) 



and the real roots, as before, as 



z 



± 



-q''c±e~ 



vt > 0, 



0, . . . ,m 



1. 



(4.9) 



The choice of branch cuts in the counting function (3.6), (3.7) is suitable for the 
forward bias regime. In the reverse bias regime we must redefine the counting function. 
The form we use is 



m—l 
I I 

iYLiz 

1=0 
L-2m-l 

1=1 



where now 



9{z) 
9h{z) 
K{w,z) 
K,{w,z) 



In 
In 

-In 
-In 



i9h{z) + iYl (^^(^r, ^) + Kr{4^ ^)) 

zil - qzy\ 

""-h + a 1 + z/c 



(4.10) 



(1-^) 
1 — q'^z 



z{l-z^ 
1 — qz/w 1 



+ 2 In c 



qm-l _|_ qo^^ I _|_ g^~; 



q^zw 



1 — qw/ z 1 — zw 
1 — q'^wjz 1 — c^wz 



1 — qwj z \ — wz 
With this definition, the complex roots satisfy 

27r, 



Inw, 



— In(g2;). 



(4.11) 



y^) 



'J' 



2 



L-2m- 1, 



with the integers L- that correspond to the first excited state given by 



The real roots satisfy 



L 



— m, + J. 



27r 



YM) = ^4^ 



k = 0, 



,m - 1, 



(4.12) 



(4.13) 



(4.14) 



but the I^ are not distinct. To find numerical solutions, we use an alternative form of 
the counting function that numbers all roots consecutively (see Appendix A). However, 
the form (4.10), (4.11) is better suited to the following asymptotic analysis. 
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4-2.1. Counting function for the complex roots We can account for the contribution 
of the isolated roots in (4.12) and obtain Bethe equations for the complex roots alone. 
Ignoring exponentially small parts, the sum over real roots in (4.10) telescopes, 

m—l 



1=0 



'-0 ^ ^ 

2m ln{qz) + 2 In 
— ln{z + c) — In 



1 + 



m+l 



qm-1^ 



(c + qz) 



l + qm+icz){l + q'^cz) 



(4.15) 
[l + qcz){l + cz) ^ ' 

The complex roots lie on a simple curve in the complex plane, which approaches 
a closed contour as L — )■ oo. From (4.2) and (4.4), we see that the size of the contour 
scales as g"*~^. We therefore change to the scaled variables 

z = g™-'C, (4.16) 

and take into account the 2m real roots "missing" from the contour by defining a new 
counting function 

L 



Y 



L-2m 



(C) 



-Ydq'^-'O 



(4.17) 



L-2m' 

With the change to the scaled variables C,, we may drop any 0{q"^) contributions. For 
example. 



g{z) = g{q-'\) = (m - 1) Ing + In ^_^^^^ 

~ (m- l)lng + lnC. (4.18) 

Collecting terms by order of 1/(L — 2m), the new counting function is given by 



lYL-2m{0 = m + 



L-2m 



UO + 



1 



L-2m-l 



L-2m 



E ^(^''0, 



(4.19) 



1=1 



with 



^(C) 


= lnC 


^b(C) 


= -ln 


n^x) 


= -ln 



ln(C) + 21n ( ^^^"j +21n(C + gc) + 21n ("l + ^ ) ,(4.20) 



qC/^\ 



l-qco/Cj 



ln(w) 



where we recall that a is defined in (|4.4|). The telescoped sum (|4.15|) contributes terms 

(4.21) 



to both g and (?b- The Bethe ansatz equations for the scaled complex roots are 

^L-2m(Ci^ 



with 



I^ 



L-2m ■" 

L-2m 
t; — +J- 



3 



L-2m-l, 



(4.22) 
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The eigenvalue, in terms of the scaled roots C,j is 

L-2m-l 
1 (rev) 



^(re.) ^ _^(rev) _ ^ _^^^^ 



(4.23) 



where 



1=1 



£(-^) = -2(1 - g)g— 1 



qc 



m— 1 



a 1 + g'"c 



(4.24) 



includes £q and the contribution from the real roots, and the scaled bare energy is 



e(C) = -(i-g)V 



2„m-l 



c 



:i - g— ic) (1 - g™C) 



(4.25) 



5. Asymptotic Analysis 



The equations (4.21 ) are logarithmic Bethe equations for the complex roots in the reverse 



bias regime. These equations have the same form as the forward bias equations (3.6) 



with the key difference that there are now L — 2m — 1 roots instead of L — 1 (recall 



that m is defined by the condition (4.3)). The method used in 24 to derive exact large 
L asymptotics for the eigenvalue £i can be applied to the reverse bias regime. But we 
require now that L — 2m is large. 




Rei; 



Figure 8. Integration contour enclosing the complex roots. 



As a simple consequence of the residue theorem we can write 

_, L-2m-l 



L-2m 



i=i 



C1+C2 



^J{0y'L-2m{C) cot Q(L - 2m)n_2^(C) ) , (5.1) 



where C = C1 + C2 is a contour enclosing all the complex roots Q (see Figure [s]) and /(C) 



is a function analytic inside C. We will use (5.1 ) to write an integro-differential equation 
for the counting function l^L_2m(C)) but first we note some properties of l^L-2m(C)- 
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The complex roots Q are connected by the contour Im [Fi„2m(C)] = 0- This 
contour extends to a point on the negative real axis, (c- We fix C, and ^*, the endpoints 
of Ci and C2, on this contour by 

On Ci, inside the contour connecting the roots, the imaginary part of l^L-2m(C) is 
positive; on C2 it is negative. 

Using (5.1) in (4.19), we obtain a non-linear integro-differential equation for 
yL-2m{C): which, using the properties just described, is written as an integral over 
the contour of roots from ^* to ^, and correction terms integrated over Ci and C2: 

iYL-2UC) = m + ;^3^^b(C) + ^ / Kiu,C)Y'L_2J^)du 

doo 



2-n J(j 1 — Q~i{L-2m.)Y L-2m{i^) 



The correction terms can be approximated as in 1124] giving 

iYL-2m{C) = m) + ;^3^^b(C) + ^ r K{ujX)Y'L-2m{^)^^ 

■^ Cc 

12(i - 2m)' l,F;_„„(e) n-^JOJ V(i - 2m)^; ' ^ ' ' 

Here, the derivatives of K are with respect to the first argument, and we have extended 
the integration contour beyond the endpoints, so that it pinches the negative real axis 
at points C^ = Cc ± iO. 

Equation (5.4) is solved as follows. We Taylor expand the integrands of the integrals 
from ^* to C,~ and from (^+ to ^. Then we substitute the expansions 

00 00 

YL-2rn{0 = ^^(L - 2m)-"y„(C), i = Cc + ^^(^ - 2m)-"(5„ + zr/„) (5.5) 

ra=0 n=l 

and expand in inverse powers of L — 2m. This yields a hierarchy of integro-differential 
equations for the functions yn(C) 

1 f^^ — 

2/n(C) = ^n(C) + W^ / K{UJ, C)y'n{^) du. (5.6) 

27ri ;^- 
The integral is along the closed contour following the locus of roots. The first few driving 
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terms QniC) are given by 

9o{o = -im, 

9iiC) = -i^b(C) + f^i + AiF(Cc, C), 

^2(0 = ^2 + A2^(Cc, C) + /i2F'(Cc, C), 

^3(0 = '^3 + A3^(Cc, C) + fisK'iCc, C) + 1^3K"{Cc, 0- 



The functions g{() and ghiC) are defined in (4.20) and (4.20), and 
arises when we evaluate K{(^, () and take the hmit (^ — )■ Cc- 



(5.7) 



(5.8) 



The coefficients k„, A„, /i„ and Un are given in terms of (5„, i/^, defined by (5.5) 



and by derivatives of j/„ evaluated at Cc- Explicit expressions are given in Appendix B 



We note that although the definition of the counting function has changed, and despite 
expanding in L — 2m instead of L, the form of these coefficients and the driving terms 



are unchanged from 24 . 



5.1. Solving for the y„ 



We may now solve (5.6) order by order. Working with the scaled roots, with the 



restrictions on parameters of the reverse bias regime we can assume that —qc lies inside 
the contour of integration, and that — -^ and — c lie outside it. The calculation follows 

° ^ qa 

in the same way as in [24], so we give here only the result for the counting function: 

1 _. 1 



YL-2miO = yoiO + 



L-2m 



yiiO 



(L 



where 



yoiO 
yiiO 



+ 

— iln 

— iln 
-2iln 
-2iln 



1 



(L - 2m) 



c 

Cc 

(_ 

Cc 



MO + o 



2m)2 
1 



1/2(0 



(L - 2m)4 



(5.9) 



(5.10) 



2iln(a) - Ai In(-Cc) + ki 



i-(lc/C,<l)o 
[-Qc/Cc, q)c 

(-C/c;g)oo 

[-Cc/c; q)oc 



+ Ailn 



(gCc/C; g)c 
(gC/C; g)c 



(5.11) 



1/2(0 
2/3(0 



K2 - A2ln(-Cc 



Cc 
Ai3 



+ A2ln 



/t3- A3ln(-Cc 

+ /i3^i(C|g)-^^3^2(C|g) 



Cc Cc' 



(gCc/C; g)c 
(?C/Cc; q)c 
(gCc/C; q) 



+ /i2^i(C|g) 



As In 



(?C/Cc; g)c 



(5.12) 
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Here (a; g)oo denotes the q-Pochhammer symbol 



g)oo = n(i 



aq 



fc=0 



and we have defined functions 

^kiC\q) = MC\q-') - MCc\q-') - iMC\q) - MCdq)) 
witlf] 

oo ^ 



n=0 



ic-q^+Kr ■ 



(5.13) 



(5.14) 



(5.15) 



5.2. Boundary conditions 



Substituting the expansions (5.5) into the boundary conditions (5.2), which fix the 



endpoints ^ and ^*, we obtain a hierarchy of conditions for yn{Cc), e.g. 



YL-2UO = yoiO + 
= yoiCc) + 



L-2m 



1/1 (0 + 



(L - 2m)2 



1/2(0 + 



L-2m 



bi(Cc) + yd(Cc)(5i + ir?i)] + ... 



n 



n 



L-2m 

Solving this equation order by order, we find 

1 



and 



Ai = 2i, Cc = -: 



A3 = Ai2 = A2 = Ki = K2 = 0, 



^^3 = Cc/^S = -ITT Cc- 



5. 3. Eigenvalue with largest non-zero real part 



(5.16) 



(5.17) 



(5.18) 



As was done for the counting function, we use (5.1) in (4.23) to obtain an integro- 



differential equation for the eigenvalue. Evaluating the resulting integrals in the same 
way as for the counting function itself, we obtain 

L-2m 



£ = -£, 



(rev) 



27r 



"^Cc n>0 



where the integral is over the closed contour on which the roots lie, £q and e are 



given in (4.24) and (4.25), and 



Co = Ai£(Cc), 

ei = A2e(Cc) +Ai2^(Cc), 

62 = A3£(Cc) + Ai3^(Cc) + 1^3^" iCc)- 
The definition of V'fc(Ck) used in i24i is incorrect for expressions such as V'fe(C^^k^^) 



(5.20) 



17 

Substituting the expansion (5.9) for YL_2m{C) iiito (5.19) we arrive at the following 



result for the eigenvalue of the transition matrix with the largest non-zero real part 

^r- ,fl.,::i}:^l ^o(,^±:^. (5.21) 



(L-2m)2g2(m-i)_^2 ' \{L-2m)^ 
Making the substitution a = q"^~^a, we see that this is 

,(rev) _ 1 7r2(l - q) 



^'^ - (L- 2m)2 a-^-a ^ ^ [{L - 2mf ' ' ^^'^^^ 



which differs from the expression in the forward bias regime (3.10) only by changing 



L — 7- L — 2m. In Section [7] we discuss the physical interpretation of this result. 
6. Crossover region 



(is) 

We have found two very different expressions for the eigenvalue: £{ in (4.7) describing 
the relaxation for relatively small L, dominated by the isolated real roots, and £| in 



(5.21) describing the asymptotic relaxation when L — 2ra is large. As L increases, the 
contour of complex roots begins to form. The first approximation is no longer valid, 
and the second cannot be applied until the contour is sufficiently large. 

Nevertheless, we can compare the magnitude of the two expressions in this region. 



Assuming m' = m, we can approximate (4.7) as 

^f)^-g'"-i(l-g)Q-(l + g)gc^. (6.i; 

The asymptotic expression (|5.21[) can be approximated as 



^r. _,»-.(! _,)_^i. (6.2) 

The contour begins to form when L — 2m = 2, and at this point we see that the two 
expressions are comparable. Figure |9] shows the two expressions and the numerically 
calculated value in the crossover region. 

7. Conclusion 

We have computed the relaxation rate on the coexistence line of the partially asymmetric 
exclusion process with reverse biased boundary rates. Such rates introduce a length scale 



m, !^ — log(— ac)/logg in the system, determined by (4.3). There are two distinct sub 
regimes. When the system size is small compared to the reverse bias, i.e. L < 2m, the 
relaxation is exponential in the system size and the inverse relaxation time is given by 



(4.7). For large systems, L 3> 2m,, the relaxation on the coexistence line is diffusive but 



the inverse relaxation time vanishes with the square of a reduced system size L — 2m. 



This behaviour can be obtained by an effective one particle diffusion model as in 23 , 35 
but in a reduced system of size L — 2m. 

We thus observe the effects of the reverse bias on the relaxation rate even with 
non-zero (but small) rates in the forward direction - provided the rates in the reverse 
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CO -4x10" 



-8x10-'" h 

20 40 60 80 

L 
Figure 9. E{^ (dashed), E{^'^^ (dotted), and numerical value (datapoints connected 



by thin solid line) for m — 20 and parameters (4.5 1 



direction are large enough. We suspect that in this case the forward contribution to the 
stationary current of particles is not strong enough to destroy the reverse bias phase, 



contrary to the argument in 13 . It would be of interest to revisit the calculation of the 
stationary current to see whether this is indeed the case. 

The corresponding physical interpretation of these results is the following. In the 
setup of this paper, where the bulk bias is from left to right, the system will fill up 
completely on the right hand side over a length 2m. In the remaining space of size 
L — 2m, the system behaves as a regular PASEP on the coexistence line, and forms 
a uniformly distributed domain wall. As far as we are aware, the stationary density 
profile of the PASEP in the reverse bias regime has not been computed analytically, and 
it would be interesting to see whether this picture is confirmed by such a calculation. 
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Appendix A. Numerical solutions of the Bethe equations 

The logarithmic form of the Bethe equations is numerically stable for large L, as long as 
an appropriate form of the counting function is used. In the forward bias regime, that 



form is (3.6). For the reverse bias regime the situation is more complicated. The form 
(4.10) is suited to the asymptotic analysis - by assuming the positions of the isolated 
roots, we choose the branch cuts of the counting function to number the complex roots 
consecutively. However, to find numerical solutions, we arrange for the counting function 
to number all roots consecutively. That is, we define a counting function 1^,^' so that 
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the Bethe equations become 






ZTT 
yL,m'{Zj) — Ij, j — 1,. . . ,L 1, 


(A.l) 


with 








/, = -L/2+j. 


(A.2) 


The complej" 


: contour is formed by the roots 






Zj, j = 1 . . . L/2 - m' - 1, L/2 + m',...,L-l, 


(A.3) 


while the isolated roots are 






Zj = z,=-q'c-e-''^\ j = L/2-k-l, 


(A.4) 


and 








h = 4 = -^'c + e-"^'', J = L/2 + k, 


(A.5) 


for A; = 0, . . 


m'-l. 





We can think of (4.10), with m = (no isolated roots), as li„m'=o, which then is 
modified as m' increases. For example, the first pair of isolated roots are z^. We work 
directly with the v^ to reduce the required floating point precision, and for Zq rewrite 
the boundary term 

In (4 + c) = In r (-C + e-^o^) + c) = -u+L. (A.6) 

For Zq we must also consider the phase as 

In {zq + c) = In (-e-"^^) = -u^L ± ivr. (A.7) 

Referring to Figure 5(b)[ Zq comes from the pushed out complex root with negative 
imaginary part, which tends to 

-c-e-^o^ -iQ-. (A.8) 

Therefore, the appropriate choice of phase is 

In [zq + c) = -UqL - ivr. (A. 9) 

The interaction terms from K{w, z) must also be handled similarly. Roots z^ and 
z^_-^ differ by a factor of g, apart from the exponentially small correction. This gives 
terms such as 

In (±e-''i^ T e^'^^-i^) . (A. 10) 

When the argument of the logarithm is negative jjj] we must choose the appropriate 
phase. Again, we consider how the isolated roots form, but must now consider their 
phase relative to all other roots. 

With this form of the counting function, we have found numerical solutions for 
large L and m' using standard root finding techniques. The key difficulty that remains 
is finding the transition point from a pair of complex roots pinching the real axis (Figure 
5(b)[ ) to a pair of isolated real roots (Figure 5(c) ). 



Note that we assume that vf > vt, for k < k' , based on the numerical solutions. 
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Appendix B. Expansion coefficients 



In this appendix we list tlie coefficients arising in the expansion (5.6), (5.7) of the integral 
equation for the counting function YL_2m{C)- ^^ the following list we abbreviate y^(Cc) 
by y'n- We note that by definition 6n and //„ are real quantities. 

«i = -Vo^u (B.l) 

Ai = y'o-, (B.2) 

71 

«2 = - y'oS2 - y'A - Ivoi^l - ril), (B.3) 

A2 = - (1/0^2 + y'lVi + 1/0^1^1) , (B.4) 

/^2 = 1/0 ' (B-5) 

TT 

/^3 = - y'oh - y'i^2 - Z/2^1 - yl{h^2 - r]ir]2) - 2^1 ('^1 " ^1) 

-^l/o%('5?-3r/D, (B.6) 

A3 = - (2/0^3 + 1/^2 + y'2'ni + Z/o ('^1^2 + 52?7i) + 5iriiy'l 

+ k'V(3<5?-^?)V (B.7) 



6 

1 



/^3 = - 

TT 



2 2 // 

.2 ^IN , ^1/0' 



y^{^\ri2 + 52??i) + y'lf^ir/i + 2/0^1(^1 ~ y ) + y^^i 



(B.8) 



^3 = ^fy;r^i(35?-r^D-vr^^). (B.9) 

67r V 2/0/ 



References 



[1] 
[2] 



[3] 
[4] 
[5] 

[6] 
[7] 

[8] 



Spitzer F 1970 Interaction of Markov processes Adv. Math. 5 246-290 

Liggett T M 1985 Interacting particle systems {Grundlehren der Mathematischen Wissenschaften 

[Fundamental Principles of Mathematical Sciences] vol 276) (New York: Springer- Verlag) ISBN 

0-387-96069-4 
Derrida B 1998 An exactly soluble non-equilibrium system: the asymmetric simple exclusion 

process Phys. Rep. 301 65-83 
Schiitz G M 2001 Exactly solvable models for many-body systems far from equilibrium Phase 

transitions and critical phenomena, Vol. 19 (San Diego, CA: Academic Press) pp 1-251 
Golinelli O and Mallick K 2006 The asymmetric simple exclusion process: an integrable 

model for non-equilibrium statistical mechanics /. Phys. A 39 12679-12705 [Preprint 

|arXiv:cond-mat/0611701 | 



Blythe R A and Evans M R 2007 Nonequilibrium steady states of matrix-product form: a solver's 

guide J. Phys. A 40 R333-R441 {Preprint arXiv: 0706. 167 8) 
Schmittmann B and Zia R 1995 Statistical Mechanics of Driven Diffusive Systems Phase 

Transitions and Critical Phenomena vol 17 ed Domb C and Lebowitz J (London: Academic 

Press) 
Derrida B, Evans M R, Hakim V and Pasquier V 1993 Exact solution of a \D asymmetric exclusion 

model using a matrix formulation J. Phys. A 26 1493 



21 

[9] Schiitz G and Doniany E 1993 Phase transitions in an exactly soluble one-dimensional exclusion 
process J. Stat. Phys. 72 277 {Preprint 'arXiv : cond-mat/9303038) 
[10] Sandow S 1994 Partially asymmetric exclusion process with open boundaries Phys. Rev. E 50 



2660-2667 [Preprint jarXiv : cond-mat/9405073 1 



[11] Essler F H L and Rittenberg V 1996 Representations of the quadratic algebra and 

partially asymmetric diffusion with open boundaries J. Phys. A 29 3375 [Preprint 

arX iv:cond- mat/9506131) 
[12] Sasamoto T 1999 One-dimensional partially asymmetric simple exclusion process with 

open boundaries: orthogonal polynomials approach J. Phys. A 32 7109 [Preprint 

arXiv : cond-mat/99 10483 ) 
[13] Blythe R A, Evans M R, Colaiori F and Essler F H L 2000 Exact solution of a partially 

asymmetric exclusion model using a deformed oscillator algebra J. Phys. A 33 2313 [Preprint 

arXiv : cond-mat/99 10242 ) 
[14] Schiitz G 1997 Exact solution of the master equation for the asymmetric exclusion process J. Stat. 



Phys. 88 427-445 [Preprint arXiv : cond-inat/9701019 1 



[15] Johansson K 2000 Shape fluctuations and random matrices Comm. Math. Phys. 209 437-476 

[Preprint arXiv:math/9903134) 
[16] Prahofer M and Spohn H 2002 Current fluctuations for the totally asymmetric simple exclusion 

process Progress in Probability vol 51 ed Sidoravicius V (Boston: Birkhauser) p 185 [Preprint 

'arXiv : cond-mat/0101200 ) 
[17] Tracy C and Widom H 2008 Integral formulas for the asymmetric simple exclusion process Comm. 



Math. Phys. 279 815-844 [Preprint [arXiv : 0704 . 2633 1 



[18] Dhar D 1987 An exactly solved model for interfacial growth Phase Transit. 9 51 

[19] Gwa L H and Spohn H 1992 Six-vertex model, roughened surfaces, and an asymmetric spin 

Hamiltonian Phys. Rev. Lett. 68 725-728 
[20] Gwa L H and Spohn H 1992 Bethe solution for the dynamical-scaling exponent of the noisy Burgers 

equation Phys. Rev. A 46 844-854 
[21] Kim D 1995 Bethe ansatz solution for crossover scaling functions of the asymmetric XXZ 

chain and the Kardar-Parisi-Zhang-type growth model Phys. Rev. E 52 3512-3524 [Preprint 

arXiv:cond-mat/9503169l 
[22] de Gier J and Essler F H L 2005 Bethe Ansatz Solution of the Asymmetric Exclusion Process with 

Open Boundaries Phys. Rev. Lett. 95 240601 [Preprint arXiv : cond-mat/050870 7) 
[23] de Gier J and Essler F H L 2006 Exact spectral gaps of the asymmetric exclusion process with 

open boundaries JSTAT 2006 P12011 [Preprint arXiv : cond-mat/0609645') 
[24] de Gier J and Essler F H L 2008 Slowest relaxation mode of the partially asymmetric exclusion 



process with open boundaries J. Phys. A 41 485002 [Preprint arXiv: 0806. 3493) 



[25] M Uchiyama, T Sasamoto and Wadati M 2004 Asymmetric simple exclusion process 

with open boundaries and Askey-Wilson polynomials J. Phys. A 37 4985 [Preprint 

arXiv : cond-ma t/03 12457 ) 
[26] Uchiyama M and Wadati M 2005 Correlation function of asymmetric simple exclusion process with 

open boundaries J. Nonlinear Math. Phys. 12 676 [Preprint arXiv : cond-mat/0410015) 
[27] Cao J, Lin H Q, Shi K J and Wang Y 2003 Exact solution of XXZ spin chain with unparallel 

boundary flelds Nucl. Phys. B 663 487 - 519 [Preprint arXiv:cond-mat/021216 3) 
[28] Nepomechie R I 2003 Functional Relations and Bethe Ansatz for the XXZ Chain J. Stat. Phys. 

111(5) 1363-1376 [Preprint arXiv:hep-th/0211001) 
[29] Yang W L and Zhang Y Z 2007 On the second reference state and complete eigenstates of the 



open XXZ chain Journal of High Energy Physics 2007 044 (Preprmt |arXiv:hep-th/0703222| 



[30] Simon D 2009 Construction of a coordinate Bethe ansatz for the asymmetric simple exclusion 
process with open boundaries Journal of Statistical Mechanics: Theory and Experiment 2009 
P07017 [Preprint [arXiv : 0903 . 4968 ) 

[31] Crampe N, Ragoucy E and Simon D 2011 Matrix Coordinate Bethe Ansatz: Applications to XXZ 



22 



and ASEP models ArXiv e-prints {Preprint arXiv: 1106.47121 
[32] Proeme A, Blythe R A and Evans M R 2011 Dynamical transition in the open-boundary 
totally asymmetric exclusion process Journal of Physics A Mathematical General 44 035003 — h 



{Preprint arXiv : 1010 . 5741 ) 



[33] Crampe N, Ragoucy E and Simon D (2010) Eigenvectors of open XXZ and ASEP models for a 
class of non-diagonal boundary conditions J. Stat. Mech. P11038 {Preprint arXiv: 1009.4119) 

[34] de Gier J and Essler F H L 2011 Large deviation function for the current in the open asymmetric 
simple exclusion process Phys. Rev. Lett. 107 010602 {Preprint arXiv: 1101 .3235| 

[35] Kolomeisky A B, Schiitz G M, Kolomcisky E B and Straley J P 1998 Phase diagram of one- 
dimensional driven lattice gases with open boundaries Journal of Physics A: Mathematical and 
General 31 6911 



